program main

    use easy_stats_m, only: cov
    
    print *, cov([real(8) :: 1,2,3,4,5])

end program main

! 2.00000000
